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Abstract 

We calculate the electromagnetic form factors of the nucleon to third chiral order in manifestly 
Lorentz-invariant effective field theory. The p and u mesons as well as the A(1232) resonance are 
included as explicit dynamical degrees of freedom. To obtain a self-consistent theory with respect 
to constraints we consider the proper relations among the couplings of the effective Lagrangian. For 
the purpose of generating a systematic power counting, the extended on-mass-shell renormalization 
scheme is applied in combination with the small-scale expansion. The results for the electric and 
magnetic Sachs form factors are analyzed in terms of experimental data and compared to previous 
findings in the framework of chiral perturbation theory. The pion-mass dependence of the form 
factors is briefly discussed. 

PACS numbers: 12.39.Fe, 13.40. Gp, 14.20.Dh 
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I. INTRODUCTION 



Electromagnetic form factors parameterize the single-nucleon matrix element of the elec- 
tromagnetic current operator and provide important information about the structure and 



composition of the nucleon (see, e.g., Refs. [l|-|3j for an overview). Furthermore, they are 
important input to high-precision tests of quantum electrodynamics as well as the Standard 
Model of particle physics. In the space-like region, the proton form factors have been mea- 
sured with great accuracy over a wide range of momentum transfer in experiments on elastic 
electron-nucleon scattering. Neutron form factors are not as well known since they have to 
be extracted from scattering experiments with deuterium or 3 He. Despite the wealth of the 
available data there are still open issues such as the value of the proton charge radius deter- 
mined from the Lamb shift in muonic hydrogen on the one hand [4J and electronic hydrogen 
Lamb shift measurements and elastic electron-proton scattering on the other hand [5]. A 
new generation of precision measurements of electromagnetic form factors at low momentum 
transfer has been and is presently performed at the Mainz Microtron MAMI @]. 

Chiral perturbation theory (ChPT) j8-lo|| is the effective field theory (EFT) of quan- 
tum chromodynamics in the low-energy domain (for an introduction and review see, e.g., 
Refs. (ill 12|). The first form factor calculation was performed in the early relativistic ap- 
proachjlO||. in which, however, the power counting of low-energy dimensions was still an 
open issue due to the additional heavy-mass scale introduced by the nucleon. Later on, 
the problem of setting up a consistent power counting in EFT with heavy degrees of free- 
dom was handled by employing the heavy-baryon approach [l3|, 14] and, more recently, by 



choosing suitable renormalization prescriptions in a manifestly Lorentz-invariant framework 
15Hl8|. Within the heavy-baryon approach, calculations of the form factors were performed 



in Refs 
3, in Ref. 



14j . 1 1911 and, including the A (1232) resonance in terms of the small-scale expansion 



Applying different renormalization schemes, form factor calculations have 



23 and 



25 



been performed within manifestly Lorentz-invariant baryon ChPT up to 0(q 4 ) 
to 0(q 3 ) including the leading-order corrections due to the A resonance j24[. In general, such 
calculations describe the experimental data only for a small range of momentum transfer 
(Q 2 < 0.1 GeV 2 ). On the other hand, the p, u, and <p mesons were included dynamically in 
the effective Lagrangians of Refs. 
which 



A systematic re-summation of higher-order terms 
in an ordinary chiral expansion, would contribute at higher orders beyond 0(q 4 ), 
results in an improved description of the data even for higher values of Q 2 , as expected on 
phenomenological grounds. The re-organization proceeds according to well-defined rules [26[ 
so that a controlled, order-by-order calculation of corrections is made possible. 

A covariant formalism for massive vector fields involves Lagrangians with constraints, 
because one typically introduces unphysical degrees of freedom 
Refs. (T 



25 



271 ] . In comparison with 



the present paper considers the conditions on the form of the Lagrangian 

We use 



28, 29 



imposed by the demand for a self-consistent theory in terms of constraints 
the extended on-mass-shell (EOMS) scheme [l8[ to generate a systematic power counting in 
the presence of heavy degrees of freedom 26[. As a result, we obtain an effective Lagrangian 
which is renormalizable in the sense of effective field theory [30| and which is consistent with 
the constraints order by order. 

Because of its close proximity to the ground state and its strong coupling to the pion- 
nucleon-photon system we also include the A(1232) resonance explicitly in our effective 
theory. The nucleon-A mass splitting is treated as an additional small parameter (small- 
scale expansion (2(|). Akin to the vector- meson case we respect the constraints on the 
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possible interaction terms to obtain a self-consistent theory describing the right number of 



degrees of freedom [3ll . 132 



This paper is organized as follows. In Sec. II the definitions of the Dirac and Pauli as 
well as the Sachs form factors are given. We briefly discuss those elements of the most 
general effective Lagrangian relevant for the subsequent calculation and state the applied 
power-counting rules in Sec. III. In Sec. IV we discuss the fit of our results to experimental 
data. The final results for the Sachs form factors are presented and analyzed. Section V 
contains a short summary. 



II. ELECTROMAGNETIC FORM FACTORS OF THE NUCLEON 

Neglecting the contributions due to heavier quarks, the electromagnetic current operator 
is given by 

J"(x) = |ii(s)7Ms) " \d{x)^d{x) = q(x) Q + |) ^q(x), (1) 
and the interaction with an external electromagnetic four-vector potential reads 

•^e.m. &J^Afj,i (2) 

where e > denotes the elementary charge. In the one-photon-exchange approximation, the 
electromagnetic form factors are defined via the matrix element 



{N(p f )\J*{0)\N(p i ))=u(p f ) 



7^(Q 2 ) + ^ F 2 N (Q 2 



u{pi), N = p,n, (3) 



where m p denotes the proton mass, q = pf — pi is the four- momentum transfer, and Q 2 = 
—q 2 > 0. The functions F^(Q 2 ) and F^(Q 2 ) are called Dirac and Pauli form factors, 
respectively. At Q 2 = 0, the Dirac form factor takes the value of the electric charge in 
units of the elementary charge and the Pauli form factor takes the value of the anomalous 
magnetic moment in units of the nuclear magneton: 

Ff (0) = 1, iff (0) = 1.793, (4) 
F™(0) = 0, F 2 n (0) = -1.913. (5) 

Our final results will be displayed in terms of the electric and magnetic Sachs form factors 



33| since these are better suited for the analysis of experimental data. The Sachs form 



factors are related to the Dirac and Pauli form factors as follows: 

g e(Q 2 ) = F^(Q 2 ) - ^F 2 N (Q% (g) 
G%(Q 2 ) = F 1 n (Q 2 ) + f/(Q 2 ). 

Sometimes it is more convenient to work with the isoscalar and isovector form factors defined 
as the sum and difference of the proton and neutron form factors, respectively, 

F^ = Ff + F?, F^ — Ff — F", i = l,2. (7) 
The isoscalar and isovector Sachs form factors are defined accordingly. 
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III. EFFECTIVE LAGRANGIAN AND POWER COUNTING 



A. Non-resonant Lagrangian 

The non-resonant part of the effective Lagrangian consists of a purely mesonic part and 
a part describing the interaction of pions and nucleons. From the mesonic sector only the 
lowest-order Lagrangian £2, including the coupling to an external electromagnetic four- 
vector potential in terms of the isovector field = —eA fl r 3 /2, is needed (9[, 

F 2 F 2 M 2 F 2 

C 2 = — Tr (dpUdPtf) + ^^Tr {U ] + U) + %— Tr [{d»Urf + d^U)v,] . (8) 

Here, F denotes the pion-decay constant in the chiral limit, F n = F [1 + 0(rh)] = 92.2 MeV, 
and M 2 = 2Bm is the squared pion mass at leading order in the quark-mass expansion. In 
the isosp in-symmetric limit rh = = ma, and B is related to the scalar singlet quark 
condensate (qq) in the chiral limit (1, The pion fields are contained in the unimodular, 
unitary matrix U : 

/ d>(x)\ 

U(x) = u 2 (x) = exp f i— — j , (f> = (j) k T k . 

Collecting the proton and nucleon fields in the isospin doublet the lowest-order ttN 
Lagrangian is given by [lOj 

= ¥ (1$ - m + ^Vt5« m ) tt, (9) 



with 



= (0, + T M - «;«) V 



u 



where = — e^4 M /2 In Eq. Q, m and g A denote the chiral limit of the physical nucleon 
mass and the axial- vector coupling constant, respectively. 

The complete Lagrangians at second and third order can be found in Ref. j35|. We only 
display those terms needed for our calculation, 



. w 2 7 „. (10) 

4v = i^ d ^ D ^ fU D "^ + H - c - + ^ f (^#) D ^ + H - c - + • • • > 

where H.c. refers to the Hermitian conjugate and 

f% = u fL^uU j ± u j f Rflu u, 
fRfiu = dfJy - d v r^ - i [r M! r v \ , 



v (s) _ a (a) _ q (s) 



with r M = 1^ = -e^4 M r 3 /2. 
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B. Lagrangian containing vector mesons 



The p-meson triplet consists of a pair of charged fields, p^ = (pi M T^P2/ i )/v / 2, and a third 
neutral field, p° = ps^. Using a covariant Lagrangian formalism, self-interacting massive 



vector fields are subject to constraints [27 



It was shown in Ref. [28[ that the requirement 
for a quantum field theory of vector mesons to be self consistent in terms of constraints 
and perturbative renormalizability leads to relations among the coupling constants of the 
Lagrangian. Eventually, at leading order the self-interacting part of the most general effective 
Lagrangian for p mesons reduces to a massive Yang-Mills structure j^. 1 



where 



Cpe« = -^Tr (p^pT) + M p 2 Tr(p^), 



Tk 

Pfj,u = d^pv - dyp^ - ig[pp, p v 



(11) 



The Lagrangian of Eq. ( ITT]) contains two parameters, namely, the p-meson mass M p (in 
the chiral limit) and a coupling strength g. Under the pair of local chiral transformations 
(Vl, Vr), we choose the p mesons to transform inhomogeneously [36| (model III of Ref. (37|), 



g 



(12) 



where 



K(v L ,v R ,u) = Jv R uvl v r Vu. 



Equation (Tl2l implies p pu t— > Kp^K^ . The mass term remains chirally invariant through the 
replacement p^ — > p M — (i/g^^, the result of which transforms homogeneously under local 
chiral transformations. Neglecting terms irrelevant for the calculation of the form factors, 
the effective chiral Lagrangian can be written as 



TTp 



-Tr (p^p 



g 



+ Y Tc Mn+-~- (1 3 ) 



Besides the proper relations among the self couplings of the p mesons, the Lagrangian of 
Eq. ( fl3l) gives rise to vector-meson dominance in the sense that, both, the pm'n and the p7 
coupling contained in the mass term [second term on the right-hand side of Eq. ( TT3]) ] are 
of leading order. The d x term parameterizes a deviation of higher order. In model III of 
Ref. [13] it is neglected, i.e., set to zero. According to Ref. 37], the Lagrangian of Eq. f[T5]) 
is obtained from the most general one by performing a field redefinition and implementing 
certain relations among the coupling constants. These relations exactly correspond to those 
derived in Ref. 28j] rendering Eq. (fT3|) consistent with the constraints and renormalizable in 
a perturbative sense. 



From SU(2)-symmetry considerations alone, the interaction Lagrangian would contain one three-vector 
and two four-vector interaction terms with, in total, three independent coupling constants (see Eq. (44) 
of Ref. 0). 
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In addition to p mesons, we also include the oj meson as a dynamical degree of freedom. 
For our calculation, from the leading-order Lagrangian 37j we only need the coupling of the 
oj meson to external fields: 



4 3 2 = -/« (d^ - sruf) «w + • • • . 



(14) 



Finally, we require the coupling of vector mesons to the nucleon which for our purposes 
is given by 



r 



VN 



9 Pn 



9 



+ 7;9uj^ix 



G n ^ 



fJtV 



(J 



(15) 



Here, we have applied the universality of the p- meson coupling g P NN — 9- I n the realization 



of Ref. 

context 



36 



the universal coupling is a consequence of chiral symmetry. In the present 
it is more likely to be a consequence of consistency conditions imposed by the 
demand of perturbative renormalizability |29|. A coupling of the u meson to the nucleon 
proportional to a^ iv is not needed at third chiral order, because there is no coupling at 
leading order in Eq. f TT4"|) as opposed to the leading-order p7 coupling of Eq. f fT3"j) . 



C. Lagrangian containing the A(1232) resonance 

will be described by a vector-spinor isovector- 



The A(1232) resonance I(J P 
isospinor with components 



3/3 + N 
2 V2 I 



I- 1 - 1 



H = 0,1,2,3, 



1,2,3. 



The physical A consists of an isospin quadruplet, whereas the description above involves six 
isospin components. In order to project onto the physical deg rees of freedom, we introduce 
the isospin projection operators (see, e.g., Sec. 4.7 of Ref. [l2| for more details) 

I 1 



ij,a/3 ~ g V iT j)al3 



where the isovector components refer to a Cartesian isospin basis. Incorporating the pro- 
ection operators explicitly, the leading-order Lagrangian in n space-time dimensions reads 

£?i=^ASr(An)e 5 ^, 



(16) 



with 



aJSTOA, n)= - {(0 - m A )(T + iA(-fD" + f 

[(n - l)A 2 + 2A + 1] 7^7" 



+ 



+ 



n-2 
m A 

(n - 2) ! 
Si 



[n(n - I) A 2 + 4(n - I) A + n] 7V 



+ f ^75<T + f (7 M ^ + uVbs + f 7^757 l 



§3, 



(17) 
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Here, the covariant derivative is given by 



V 



where we parameterized T M = T^tu- In Eq. ( !T6|) . A ^ — ~ denotes an arbitrary real 
parameter and m& refers to the leading-order mass of the A. 

Since the Lagrangian of Eq. (fl6l) describes a system with constraints, similarly to the 
previously discussed case of vector mesons, the requirement for a self-consistent theory leads 



to relations among the coupling constants 31 

g 2 = A Ei, g 3 = 



1 + 2A + A 2 (n- r 



n 



gi 



The lowest-order 7riVA interaction Lagrangian reads [32 



£■11} A — 



u„.-w + h.c. 



(18) 



with the parameterization = u^^Tk, and g being a coupling constant. 2 Since physical 
quantities cannot depend on A [32|, we choose A = — 1 in the following calculations. 



D. Power counting 

We assign a low-energy order D to each renormalized diagram. The value of D is de- 
termined with the following power-counting rules: A pion propagator counts as 0(q~ 2 ), a 
nucleon propagator as (9(g _1 ), vertices derived from £2 count as 0(q 2 ), and vertices from 
C^ N count as 0{q l ). Both p-meson and w-meson propagators count as O(q ) while vertices 
from £l 3 2 count as Otf). From the listed terms of C np and C nVN vertices of 0{<f>) to Otf) 
can be derived. The A propagator counts as 0(q~ 1 ) and vertices from C^na coun t as 0{q). 
Finally, we assign the order 0(q) to the mass difference 5 = m^—m. In order to renormalize 
the loop diagrams in such a way that they respect the above power counting, we apply the 



EOMS scheme 18 



IV. RESULTS AND DISCUSSION 



All Feynman graphs contributing to the calculation of the electromagnetic form factors 
up to and including 0(q 3 ) are displayed in Fig. [TJ 3 The 16 loop diagrams of Fig. [TJ are 
grouped into three, independently current-conserving subsets. In the following, we refer to 
the diagrams (7)— (11) as set 1, the diagrams (12)-(18) as set 2, and the diagrams (19)-(22) 
as set 3, respectively. Set 1 consists of diagrams proportional to g^/i 7 " 2 containing only pion 



2 The sign convention in Eq. (II 8[) is chosen such that SU(4) symmetry imp lies the relations g 1 = and 
g = fv^g^ among the coupling constants of Eqs. ©, (HZ]), and (TH]) [20|. 

3 After renormalization, the diagrams (14), (17), and (18) are at least of 0(q 5 ) and can thus be neglected 
in the numerical analysis. 
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m N 


M n 


M p 


M w 






Sa 


9 


g 


0.938 


0.140 


0.775 


0.783 


1.21 


0.0922 


1.27 


5.93 


1.13 



TABLE I: Input parameters: The masses and F n are given in units of GeV; the coupling constants 
gA, g, and g are dimensionless. 

loops while set 2 consists of diagrams proportional to g 2 and g 2 containing pion loops as 
well as vector-meson loops. Finally, set 3 contains all pion-loop diagrams involving the A 
resonance, thus being proportional to g 2 . 

Summing up all contributions and multiplying them with the wave function renormaliza- 
tion constant yields the final expressions for the form factors. To render the results for 
the unrenormalized form factors finite, we apply the modified minimal subtraction scheme 
of ChPT (MS) [9(| . Beyond that, we perform finite subtractions according to the EOMS 
scheme jl8| such that the power counting of Sec. IIII Dl is respected. To the given order, 
the product of the wave function renormalization constant and the tree-order diagrams sub- 
tracts all power-counting- violating terms of the loop diagrams in the Dirac form factor F\. 
In agreement with the Ward identity, we obtain Ff(0) = 1 and F"(0) = for the proton and 
neutron, respectively. On the other hand, the loop contributions to the Pauli form factor 
F 2 contain power-counting-violating terms. All subtraction terms are analytic in the pion 
mass and momenta and can be absorbed in the renormalization of the available coupling 
constants. 



A. Fixing of the LECs 



To evaluate the form factors numerically, the parameters of the effective Lagrangian need 
to be fixed. The masses, the axial- vector coupling constant, and the pion-decay constant are 
expressed in terms of their physical values, because the difference to the respective values 
in the chiral limit is beyond the accuracy of our calculation. 

Using the KSRF relation [H, S3, 



Mi 



2g 2 F\ 



(19) 



generated by the combination of chiral symmetry and the consistency of the EFT with 
respect to renormalizability [29l ]. we obtain g = 5.93. Moreover, we take g = 1.13 as 



32 . The numerical values of the above 



obtained from a fit to the A — > ttN decay width 
parameters are summarized in Table [H 

To determine the renormalized low-energy constants Cq and C7, we fix the Pauli form 
factors F 2 (Q 2 ) at Q 2 = in accordance with Eq. (j^D- The expansions of the anomalous 
magnetic moments of the proton and neutron read 



1 ^~*p 
k„ = c-jm + 2m cq - 

n n = c 7 m - 2m ( c 6 - ^ 



8F 2 



M + 



71 



gA m 



8F 2 



M — 



71 



9F 2 7T 2 

g 2 m 
9FV 



6 



In 



M 
26 



hi" 

25 



V5 2 - 


M 2 


5 


Vs 2 - 


M 2 



ln(X) 



lnpf) 



+ ■ ■ 



(20) 



S 





(*) 

K 7TN 


0) 

K 7TN 


w 

K 7rA 


(«) 

K 7rA 


Expanded 





1.98 





-0.222 


Full 


0.169 


1.35 


-0.0120 


-0.150 



TABLE II: irN and 7rA loop contributions to the isoscalar and isovector anomalous magnetic 
moments. 



with 

X = ^ , 6 = m A -m. (21) 

In Eq. (|20|) . the ellipses refer to terms scaling at least as t 2 under M i— > tM and 5 h- > 
t<5. Because of the chosen renormalization scheme, the Feynman diagrams of set 2 do not 
contribute to the magnetic moments. The non-analytic terms for the magnetic moments of 



Eq. (120]) coincide with those of Ref. [211 ] . Using the values of Tablefl]for the input parameters, 
the irN and 7rA loop contributions to the isoscalar and isovector magnetic moments are 
shown in Table [III Keeping only the leading-order terms, the pion loop contributions are 
purely isovector as in ordinary ChPT at 0(q 3 ). On the other hand, evaluating the full 
expressions modifies the isovector pieces and also generates isoscalar contributions. 

Adjusting the complete results for the magnetic moments to their empiric values yields 

c 6 = 1.39 GeV" 1 , c 7 = -0.148 GeV" 1 , (22) 

with c 6 = c 6 — G p /(2g). In Sec. IIVCI we compare our results with explicit A contributions 
to those without the A resonance. Therefore, we also state the values for the couplings of 
the latter case, namely, 

c 6 = 1.35 GeV" 1 , c 7 = -0.154 GeV" 1 . (23) 

The remaining six free low-energy coupling constants G p , f w , g u , <i 6 , d 7 , and d x are 
determined by simultaneous fits of all four Sachs form factors to experimental data for 
different regions of momentum transfer. As the data basis for the fits, we use the extensive 
proton cross section data set from Refs. 0, S3] and the neutron form factor data from Refs. 



[4lH59|. Because of the small numerical contributions originating from the diagrams of set 
2, the fits depend only marginally on the individual values of g u and f u . On the other hand, 
the product of g w and f w stemming from the tree diagram (6) of Fig. [T]is much much more 
influential for the final result. Thus we fix f u at 0.1 and only use the product f u ■ g w as an 
independent fit parameter. The results for the fitted renormalized couplings are shown in 
Table ED 

As indicated by the respective values of the reduced chi-square test (Xrcd = X 2 /d-o.f.), 
the adjusted results including only vector mesons as explicit resonant degrees of freedom 
show better agreement with experimental data than the results incorporating also the A 
resonance. As the range of momentum transfer increases, the respective values for x^ ed of 
the fits including the A increase faster than those without the A. We will come back to this 
feature in Sees. IIVOI and HVDl 



9 



Resonances 


O 2 


9uj 


fu) ' 9ui 


G p 


d 6 


d 7 


d x 


^red 




0.2 


-1.06 


-0.106 


-4.84 


1.67 


-0.282 


-0.506 


1.50 




0.3 


— 1.27 


—0.127 


—4.05 


1.55 


— 0.233 


—0.512 


A 01 

4.21 


A, p, uj 


















0.4 


-1.92 


-0.192 


-1.98 


1.50 


-0.211 


-0.547 


25.30 


p, UJ 


0.2 
0.3 
0.4 


5.13 
4.91 
4.49 


0.513 
0.491 
0.449 


-16.90 
-17.13 
-16.58 


0.629 
0.507 
0.490 


0.0909 
0.0991 
0.0934 


-0.134 
-0.118 
-0.130 


1.45 
1.74 
3.57 



TABLE III: Comparison of the renormalized coupling constants obtained from fits of the results 
to different ranges of momentum transfer. The second set includes the A resonance as an explicit 
dynamical degree of freedom. G p is given in units of GeV -1 , d§ and dj in units of GeV -2 , and 
Q^ ax in units of GeV 2 ; the remaining coupling constants d x , g^, and f u are dimensionless. 

B. Charge and magnetic radii 

Expanding the Dirac and Pauli form factors for small values of Q 2 , 

F t (s ' V \Q 2 ) = if^(O) (l - \{{r^) 2 )Q 2 + •••), 

gives access to the mean square radii. At 0(q 3 ), the expanded mean square radii are given 
by 



((rP) 2 ) 



+ 



((rf ) ) 2 ) = -12rf 6 + 6 



M 2 32F 2 vr 2 288F 2 tt 2 
3g 2 f(M p )-g 2 J(M„) + ---, 
1 - d x g 



17 + 40 In ( - 



16F 2 7T 2 
^2 



+ 



g 



54F 2 vr 2 16 



M 2 



Si 



16F 2 7T 2 

3m 2 



16M 2 



305 In (X)! 2r/1 ,v ar/wx 
7 ==\+9 2 f(M P )-9 2 J(M^ + . 



((r?) 2 ) = + • • • 



g 2 A m 



g 2 m 



8F 2 hMk( v ) 9F 2 tt 2 VS 2 - M 2 d v ) 



In (X) + 



(7) 

60 In ( — ) + 7 





41" 




+ Y 


w 





(24) 



30 In ( — 

H J 



(25) 
(26) 
(27) 



where X is defined in Eq. ( l2"Tj) . In the case of the mean square Dirac radii, the ellipses 
refer to terms that scale at least linearly in t under M 1— > tM and 5 t— > t8. On the other 
hand, for the mean square Pauli radii, the ellipses represent terms that remain constant or 
scale with higher powers. As expected (60[, {{r^) 2 ) diverges logarithmically in the chiral 
limit, whereas {(r^) 2 } shows a 1/M singularity. The respective contributions of the vector 
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mesons are given in terms of the function 



v 



12m 4 - 37M 2 m 2 + 10M 4 (4m 4 - QM v m 2 + 5M£) In 



+ 



64m 4 (4m 2 — My) tx 2 32m 6 vr 2 

M v (-36m 6 + 70M 2 m 4 - 36M 4 m 2 + 5M 6 ) arccos (f^) 



M v f' 2 vr 2 



(28) 



(29) 



Using the couplings from the fitting procedure (see Table IITT1) . we are in the position to 
determine the numerical values for the mean square charge and magnetic radii, defined as 



32m 6 (4m 2 

which vanishes in the limit of infinitely heavy vector-meson masses, 

lim / (M v ) = 0. 

My->00 



((r P E f) 



6 dG p E {Q 2 



G P E (0) 



dQ 2 
dG p M (Q 2 



Q 2 =0 



dQ 2 



({r E f) 
((r n M ) 2 ) 



-6 



dG E (Q 2 



dQ 2 
6 dG n M (Q< 



G n M (0) dQ 2 



(30) 



Q 2 =0 



The respective empirical values are shown together with our results in Table IIVI 4 As a 
general trend we find that the proton radii are better described in terms of the calculation 
including the A resonance. In contrast, the neutron radii are in better agreement with 
the experimental results in the theory without the A resonance. In all cases, {{r 7 ^) 2 ) is 
smaller and |((^^) 2 )| is larger than their respective empirical values. The situation improves 
in the calculation featuring just vector mesons as ((r^) 2 ) is approximately in agreement 
with experiment. Again, the discrepancy especially for ((r^-) 2 ) grows towards larger Qmax- 
Even though the value of |((^^) 2 )| is notably smaller than in the results with an explicit A 
resonance, it is still larger than empirically predicted. 

In principle, one could adjust d§ and d-j to the electric radii and two LECs originating 
from the vector-meson Lagrangian to the magnetic radii, respectively. However, such an ap- 
proach would be against the purpose/spirit of introducing vector mesons, namely, generating 
curvature to extend the description to intermediate values of Q 2 . 

Finally, in Table [V] we display the individual ttN, it A, and vector-meson loop contribu- 
tions to the mean square radii. The parameters have been taken from Table [B Note that 
the o;-meson loop contributes only to the mean square radii of the proton. The contribution 
is given by the second term in the respective sum and has been evaluated with the largest 
value of Table ITTTj i.e. g w = 5.13. The difference between the full and expanded loop results 
may be regarded as an, admittedly, rough error estimate. 



C. Graphical representation of the form factor results 

Using the LECs represented in Tables HI and IHH the final results for the Sachs form factors 
are displayed in Fig. [2J Since our proton results are fitted directly to cross sections we show 



When evaluating numerical expressions, we make use of the renormalization scale [i = 1 GeV. 
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Resonances 


Qm&x 


<(^) 2 > 


({r n E ?) 




<W/) 2 > 




0.2 


0.740 


-0.288 


0.631 


0.718 


A, p, u 


0.3 
0.4 


0.744 
0.761 


-0.355 
-0.440 


0.614 
0.574 


0.700 
0.667 




0.2 


0.733 


-0.198 


0.677 


0.725 




0.3 


0.730 


-0.221 


0.672 


0.728 


p, UJ 


0.4 


0.734 


-0.252 


0.659 


0.726 


Empirical values 


0.770 


-0.116 


0.604 


0.743 



TABLE IV: Comparison of the mean square charge and magnetic radii of the nucleon obtained 
from the form factor results including and excluding the A resonance fitted to different ranges of 
momentum transfer. The mean square radii are given in units of fm 2 . The empirical values are 
taken from Ref. 51. 



Mean square radii 


((r P E ) 2 ) {{r%?) ((r p M ) 2 ) <(r^) 2 ) 


ttN, expanded 
vriV, full 


0.304 -0.0938 0.301 0.381 
0.365 -0.158 0.162 0.201 


it A, expanded 
ttA, full 


0.0312 -0.0713 0.0796 0.0517 
-0.115 0.163 0.154 0.0323 


Vector-meson loops, expanded 
Vector-meson loops, full 


0.0111 0.00554 -0.000953 0.00278 
+0.00821 +0.00311 

0.0422 -0.00897 -0.0153 0.00400 
+0.00821 +0.00311 



TABLE V: ttN, it A, and vector- meson loop contributions to the mean square radii in units of fm 2 . 
The omega loop contributes only to the proton radii and is given by the second term in the sum. 

a grey band corresponding to a direct least-squares model fit for G P E and G P M to the measured 
cross sections and thus representing the experimental data |7|]. For the neutron form factors 
our results are plotted together with the respective set of data to which they have been 
fitted to. All curves describe the corresponding experimental data reasonably well for all 
four form factors. As indicated by the values of x 2 e d m Table 1111} the description of the form 
factors adjusted to a wider range of momentum transfer worsens. While G P E and G 1 ^ are 
still in good agreement with the corresponding data for = 0.4 GeV 2 (dotted curves), 
G P M does not show sufficient curvature. The electric form factor of the neutron, G E , cannot 
be described well beyond Q 2 ~ 0.3 GeV 2 . In all cases, the slope of G E at small values of 
the momentum transfer turns out too big. 

In order to investigate the total results for the Sachs form factors in more detail, Fig. [3] 
shows the individual contributions of the different diagram sets for the fit up to Q^ax = 0-4 
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GeV 2 . The dotted lines represent the diagrams involving 7rA loops which seem to adulterate 
the results for Q 2 > 0.3 GeV 2 . Similarly to the short-dashed lines involving ttN loops only, 
the A contributions do not produce significant curvature. At larger values of Q 2 , the strong 
linear Q 2 dependence of the dotted lines, especially for G P E , G^, and G P M , together with their 
lack of curvature at 0(q 3 ) cannot be compensated by the strongly curved tree contributions 
of the vector mesons (long-dashed lines). For this reason, the quantitative description of 
the data worsens for larger values of Q 2 at the considered chiral order if the A is included 
explicitly. In agreement with Ref . [25[ , the numerical contributions resulting from the vector- 
meson loop diagrams, denoted by the dash-dotted lines, turn out to be small. 

Figure H] displays the electromagnetic Sachs form factors without explicit A resonance 
(see Table HIT)) . As opposed to Fig. [2j all fits describe the related data remarkably well for all 
four Sachs form factors. Consequently, the consistent inclusion of p and u mesons provides a 
satisfactory description of the electromagnetic form factors in the momentum transfer region 
< Q 2 < 0.4 GeV 2 already at third chiral order. In our calculation, the additional inclusion 
of the <p meson results in only a small numerical adjustment of the parameters in Table II III 
and does not generate a visible improvement of the form factor curves. 



In Ref. [25|], the electromagnetic form factors of the nucleon were calculated in Lorentz- 
invariant ChPT using the EOMS renormalization scheme up to and including 0(q A ). The 
vector mesons p, u, and were incorporated explicitly using parameterization II of Ref. (37| . 
The self-consistency relations for the p meson couplings, as discussed in Sec. Ill B, were not 
considered in the effective Lagrangian. In order to discuss to what extent the self-consistent 
inclusion of vector mesons influences the description of the form factors, in Fig. Owe dis play 
a comparison of our calculation, including explicit vector mesons only, and that of Ref. [25] . 
Even though our calculation involves one vector-meson degree of freedom less, namely the 
(f) meson, and is only up to 0(q s ), the results for all Sachs form factors are slightly closer 
to experimental data for Q 2 > 0.2 GeV 2 . In particular, our curve for the electric neutron 
form factor shows a better trend for larger values of the momentum transfer. The improved 
description can be explained by the following two reasons. On the one hand, the consistency 
relations lead to a parameterization which features the cry and pirn couplings at leading 
order in distinction to the parameterization used in Ref. [25J featuring such couplings only 
at next-to-leading order. The re-ordering of terms changes the results favorably for larger 
values of Q 2 . On the other hand, we fit our free coupling constants to the global trend of the 
form factor curves instead of adjusting them to the electric and magnetic radii and taking 
the remaining ones from results based on dispersion relations. This approach allows for a 
better overall descriptions of the form factor curves. 



D. Estimate of higher-order effects 

In order to estimate the uncertainties originating in the truncation of the expansion 
at 0(q 3 ), we add polynomials in Q 2 to the Sachs form factors. The explicit form of the 
corresponding polynomial is motivated by an analysis of the maximal powers of Q 2 in the 
isoscalar and isovector Dirac and Pauli form factors, taking into account the available chiral 
structures at 0(q 3 ) with adjustable LECs. We find that the neglected structures which 
would appear in an 0{q i ) calculation can be parameterized as 

AGf = (a N Q 2 - b N ^j Q\ AC?£ = (a N Q 2 + b N ) Q 2 , (31) 
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((r P E ) 2 ) 


({r%?) 




(W f ) 2 > 


0.761 


-0.440 


0.574 


0.667 


0.718 


-0.142 


0.708 


0.758 



TABLE VI: Mean square charge and magnetic radii of the nucleon obtained from the form factor 
results including the A resonance. The first row corresponds to the 0(q s ) results, whereas the 
second row also includes additional higher-order terms according to Eq. (|3ip . The mean square 
radii are given in units of fm 2 . 

where and (N = p,n) denote unknown coefficients. The corrections to the electric 
Sachs form factors are purely of order (Q 2 ) 2 while those to the magnetic Sachs form factors 
already start at order Q 2 . To investigate the influence of these higher-order contributions 
we add them to the original expressions at 0(q 3 ) and perform a new simultaneous fit of all 
four "improved" Sachs form factors to the data. The results of this procedure are shown 
in Fig. O Clearly, the inclusion of higher-order terms according to Eq. (n?T|) changes the 
slopes of the Sachs form factors. The reason for this can be understood if one resolves 
the physical Sachs form factors into their respective isoscalar and isovector parts. Figure [7] 
indicates that an 0(q 3 ) calculation does not generate enough curvature for both, the electric 
and magnetic, isoscalar form factors. To fit to the data over a wider range of momentum 
transfer in terms of minimizing the total x 2 ed function, the respective slopes counterbalance 
the missing curvature. This overcompensation leads to a lower value for {{r p M ) 2 ) and ((r^) 2 ), 
as well as to a higher value for \({r^) 2 ) \ than in the case with additional higher-order terms. 
A comparison between the radii obtained within the two approaches is shown in Table IVIl 



E. Pion mass dependence of the form factors 

Nucleon electromagnetic form factors have been calculated in the framework of lattice 



QCD (for a recent overview, see Refs. [6lH65| and references therein). Several collaborations 
have reached pion masses down to 270 MeV. Systematic extrapolations to the physical value 
of the pion mass and the Q 2 dependence are necessary to compare the lattice form factor 
results to experimental data. Given the manifest Lorentz covariance of our results, they 
may provide useful guidance for extrapolations of lattice simulations. In lattice calculations 
it is more convenient to work with isoscalar and isovector form factors. Simulations of the 
isoscalar form factors are numerically expensive since they involve the evaluation of discon- 
nected quark loops. For the isovector form factors the disconnected quark loop contributions 
cancel. 

In Fig. [8] we display the pion- mass dependence of the Sachs form factors in the isovector 
and isoscalar channels. The quantities are derived from the results including p and uj mesons 
as explicit resonant degrees of freedom. We stress that, at the given order of our calculation, 



the vector meson masses are independent of the pion mass (see, e.g., Refs. |66H68j for a 
discussion of the quark- mass dependence of the p- meson mass). Extrapolations of lattice 
data of Ge have been performed in Refs. [69|, 70[. The results qualitatively agree with 



ours as the values of the form factor fall off more slowly with increasing pion mass. In Ref. 
fn\ the isovector magnetic moment is extrapolated to small pion masses using results from 
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heavy-baryon ChPT [72] . The pion-mass dependence and the value in the chiral limit are 
in good agreement with our corresponding findings plotted in Fig. |9j 

In the present paper we fit the unknown LECs to experimental data at the physical pion 
mass. Alternatively, they can be fitted to lattice simulation data at different values of the 
pion mass resulting in a complete theoretical prediction of the observables. Whether this is 
a useful approach largely depends upon the range of pion masses in which the low-energy 
EFT is still applicable. The region of applicability depends on the calculational scheme and 
has yet to be studied more thoroughly. 



V. SUMMARY 

We have calculated the electromagnetic form factors of the nucleon at 0(q 3 ) in mani- 
festly Lorentz-invariant baryon chiral perturbation theory including p and u mesons as well 
as the A resonance. In terms of constraints and perturbative renormalizability, we have 
incorporated the resonant degrees of freedom self-consistently into the EFT. To generate 
a systematic power counting we have applied the extended on-mass-shell renormalization 
scheme. Two of the undetermined low-energy coupling constants have been adjusted to the 
anomalous magnetic moments while the remaining six LECs have been fitted simultaneously 
to the experimental data for different ranges of momentum transfer up to Q 2 = 0.4 GeV 2 . 
We found that the results incorporating vector mesons agree well with experimental data in 
a momentum transfer region < Q 2 < 0.4 GeV 2 while those also including the A describe 
the form factors only up to Q 2 ~ 0.3 GeV 2 . For larger values of Q 2 , notably G% and G P M 
disagree with the data. This is because the 7rA loop contributions at 0(q 3 ) feature a strong 
linear Q 2 dependence without sufficient curvature. Our results including p and u mesons at 
third chiral order agree at least as well with experimental data as the previously performed 



calculations of Refs. [22j, |23] at 0(q ). This improvement is a hint towards the importance 
of a proper consideration of self-consistency relations among the couplings of the effective 
Lagrangian. The resulting form of the vector-mesonic Lagrangian gives rise to vector-meson 
dominance at leading order and deviations thereof are pushed to higher orders. 
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FIG. 1: Feynman diagrams including vector mesons and the A resonance contributing to the 
electromagnetic form factors of the nucleon up to and including 0(q 3 ). Solid, dashed, wiggly, 
double-wiggly, and double-solid lines refer to nucleons, photons, pions, vector mesons, and the 
A resonance, respectively. The numbers in the interaction blobs denote the chiral order of the 
respective vertex. 
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FIG. 2: Sachs form factors of the nucleoli at 0(q 3 ) including p, ui, and A, fitted to different ranges 
of momentum transfer Q 2 . Full lines correspond to a fit up to Q^ax = 0-2 GeV 2 , dashed lines 
up to Qmax = 0-3 GeV 2 , and dotted lines up to Q 2 nax = 0.4 GeV 2 , respectively. The grey bands 
represent empirical fits of the form factors to the measured cross sections. 
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FIG. 3: Decomposition of the Sachs form factors of the nucleon at 0(q 3 ) including vector mesons 
and the A resonance. Full lines: total results; long-dashed lines: tree contribution; short-dashed 
lines: contribution of set 1; dash-dotted lines: contribution of set 2; dotted lines: contribution of 
set 3 of Fig. [Q respectively. 
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FIG. 4: Sachs form factors of the nucleon at 0(q 3 ) including p and w fitted to different ranges of 
momentum transfer Q 2 . Full lines correspond to a fit up to Q^ax = 0-2 GeV 2 , dashed lines up to 
Qmax = 0-3 GeV 2 , and dotted lines up to = 0.4 GeV 2 , respectively. The grey bands represent 
empirical fits of the form factors to the measured cross sections. 
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FIG. 5: Sachs form factors including vector mesons. The full lines refer to the Q(q 3 ) results of this 
work including p and u mesons and the dashed lines to the 0(q ) results of Ref. 25] including p, u, 
and 4>. The grey bands represent empirical fits of the form factors to the measured cross sections. 
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FIG. 6: Sachs form factors including vector mesons and the A resonance. The full lines refer to 
the 0(q 3 ) results, whereas the dashed lines are supplemented by additional higher-order terms 
according to Eq. (|3ip . The grey bands represent empirical fits of the form factors to the measured 
cross sections. 
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FIG. 7: Isocalar and isovector Sachs form factors including vector mesons and the A resonance. 
The full lines refer to the 0(q 3 ) results, whereas the dashed lines are supplemented by additional 
higher-order terms according to Eq. (|3ip . 
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FIG. 8: Electric and magnetic Sachs form factors in the isovector and isoscalar channel including p 
and oj for different values of the pion mass M w . The full lines refer to M n = 140 MeV, the dashed 
lines to M w = 200 MeV, the dash-dotted lines to = 300 MeV, and the dotted lines to = 400 
MeV, respectively. 
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FIG. 9: Pion-mass dependence of the anomalous magnetic moment in the isovector and isoscalar 
channels and of the mean square isovector Dirac and Pauli radii; the vertical dashed line indicates 
the physical pion mass. 
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